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ABSTRACT 

The repeated discovery of supermassive black holes (SMBHs) at the centers of galactic bulges, and 
the discovery of relations between the SMBH mass (M.) and the properties of these bulges, has 
been fundamental in directing our understanding of both galaxy and SMBH formation and evolution. 
However, there are still many underlying questions surrounding the SMBH - galaxy relations. For 
example, are the scaling relations linear and constant throughout cosmic history, and do all SMBHs lie 
on the scaling relations? These fundamental questions can only be answered by further high quality 
direct M. estimates from a wide range in redshift, before further refinements to galaxy evolution 
models can be made. In this paper we determine the observational requirements necessary to directly 
determine SMBH masses, across cosmological distances, using current M, modeling techniques. We 
also discuss the SMBH detection abilities of future facilities. We find that if different M, modeling 
techniques, using different spectral features, can be shown to be consistent, then both 30 m ground- 
and 16 m space-based telescopes will theoretically be able to sample M, ~ 10 9 M Q across ~ 95% 
of cosmic history. In addition, SMBHs as small as 10 6 M Q will be sampled at a distance of the 
Coma cluster, and SMBHs as small as 10 4 M Q will be sampled in the Local Group. However, we find 
that the abilities of ground-based telescopes critically depend on future advancements in adaptive 
optics systems; more limited AO systems will result in limited effective spatial resolutions, i.e., SMBH 
detection efficiency, and forces observations towards the near-infrared where spectral features are 
weaker and more susceptible to sky features. Ground-based AO systems will always be constrained by 
relatively bright sky backgrounds and atmospheric transmission. The latter forces the use of multiple 
spectral features and dramatically impacts the SMBH detection efficiency. The most efficient way to 
advance our database of direct SMBH masses is therefore through the use of a large (16 m) space-based 
UVOIR telescope. 

Subject headings: Astronomical Instrumentation. Astronomical Techniques. 



1. INTRODUCTION 

Direct mass estimates of supermassive black holes 
(SMBHs) are made by spatially resolving the gravita- 
tional influence of the SMBH itself. There is a constantly 
growing da tabase of such direct SMBH mass ( M,) esti- 
mates (e.g. [Graham 200c| iGultekin etllll2009ft that has 
been repeatedly used to demonstrate intimate links be- 
tween M, and rudimentary properties of the surround- 
ing host galaxy. For ex ample, M, is seen to scale wit h 
the bulge luminosity (jKormendv fc Richstond fl 995). 
the s t ellar velocity di s persio n, a* (|Ferrarese fc Merrittl 
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2000), the bulge concentration 



2001), the stellar light deficit 
120091). and potentially the dark 



matter halo mass (|Ferraresd l2002t iBaes et all [2003; 
IPizzella et all[2005l) . 

These SMBH scaling relations have had a fundamen- 
tal impact on our understanding of galaxy formation and 
evolution. For example, it is possible to use scaling rela- 
tions to non-directly estimate and extrapolate M, from 
galactic bulges in which the central SMBH's gravitational 
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influence is not resolved, i.e., distant and/or small bulges. 
Therefore, scaling relations could be used to estimate 
M, across cosmic history and constrain the black hole 
mass function (BHMF). Consequently, scaling relations 
have foste red a wealth of exciting new theoretical inves- 
tigations (ICiotti fc van Albada|[200lT : lAdams et alj|2003t 
ICattaneo et al.ll2005| ) and have pot entially placed impor- 
tant limits to evolutio nary models ([Heckman et al.1120041 : 
iWvithe fc Lo^l2005h . 

Due to its small scatter the M, — a * relation has re- 
ceived the most att e ntion . However, iMarconi fc Hunt] 
( 2003) and iGrahaml (|2007l ) have shown that with care- 
ful morphological and multi- wavelength analyses, scaling 
relation scatters can be comparable. The M, — <r* rela- 
tion, characterized as logAf. = a + /Jlogtr*, has had 
many a ttempts to fit its zero-point (a), slope (0) and 
scatter dFerrarese fc Merrittl [2000l IGebhardt et alJl2000t 
iMerritt fc" Ferraresdl2001t iTremaine et al.ll2002f ). Unfor- 
tunately, these individual fits have produced consider- 
ably different results in which (3 ranges from 4.0 to 4.9, for 
example. In addition, both the intermediate mass black 
hole (IMBH, M. < 10 5 M©) and hyper-massive black 
hole (HMBH, M. > 10 10 M Q ) regimes of the BHMF are 
almost entirely unexplored. This leaves the linearity of 
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the M % — a* relation unclear. Consequently, we do not 
know the upper and lower limits to the BHMF, or even 
if such limits exist. 

If we assume that all SMBHs lie on the M, — rr, 
relation, that the relation is linear, and that fainter 
less massive galaxies contain SMBHs rathe r than com- 
pact stellar nuclei (|Ferrarese et al.l l2006a|) . then tak- 
ing ct* = 2 0km s" 1 , a value f ound in large glob- 
ular clusters ( Mevlan et al.l 119951 ) and where IMBHs 
may reside, the (3 — 4.0 slope predicts logAf. = 
4.1(+0.4, -O.4)M while the j3 = 4.9 slope predicts 
logM. = 3.4(+0.4, -O.5)Af . Taking = 444km s" 1 , a 
value seen in brightest cluster galaxies (jSalviander et al.1 
120081) and where HMBHs may reside, the fi = 4.0 slope 
predicts logM. = 9.5(+0.2, -0.1)M Q while the /3 = 4.9 
slope predicts logM. = 9.9(+0.2, -0.2)M o . Therefore, 
high accuracy direct estimates of M. from IMBHs and 
HMBHs are of critical importance. A poorly determined 
Mm — er* relation introduces large uncertainties in extrap- 
olated values of M., and influences our understanding of 
both SMBH and galaxy formation and evolution. 

However, it must be noted that scaling relations estab- 
lished in the l ocal universe ma y have experienced cos- 
mic evolution (Tre u et "al"1l2007l and references therein). 
Therefore, if evolutionary models are to be properly con- 
strained, it must be determined whether the scaling re- 
lationships themselves evolve. Consequently, an accu- 
rate and complete sample of SMBH masses must be con- 
structed from as wide a range of redshifts as possible. 
Such a database would then shift the importance of the 
scaling relations, as extrapolation of the scaling relations 
would no longer be required to constrain evolutionary 
models. 

The SMBHs that power QSOs frequently have their 
masses es timated using metho ds, such as reverberation 
mapping (jPeterson et al.ll2004h . that are calibrated from 
the scaling relations based on the unknown geometry of 
the inner broad line region (BLR). The enormous lumi- 
nosities of QSOs make them ideal candidates for estab- 
lishing SMBH masses at high redshift. However, without 
well understood scaling relations the reverberation map- 
ping calibration remains uncertain. Direct estimates if 
Mm in QSOs, while challenging, will determine if QSOs 
follow the same scaling relations as established in qui- 
escent galaxies, and will provide an absolute calibra- 
tion for continued reverberation mapping of high red- 
shift SMBHs. In the closest (un-obscured) QSO, 3C 273, 
assuming M. = 10 9 8 M Q (jPaltani fc Turlerl |2005[ ) at 
a distance of 640 Mpc, a minimum resolution of 0'.'05 
(~ 150pc) is needed to directly determine M.. However, 
using current techniques it is also a challenge to deter- 
mine an accurate value of <7», as the QSO can drown out 
the signal from the host galaxy and fill in the absorption 
features needed to determine cr». 

The future path of SMBH investigations is likely to be 
determined by the questions posed above. These can be 
summarized as follows: 

• Do HMBHs follow the same scaling relations as de- 
fined by SMBHs? 

• Do IMBHs or compact stellar nuclei exist in fainter 
less massive galaxies? 

• Are SMBH scaling relations linear? 



• Have SMBH scaling relations evolved through cos- 
mic history? 

• What are the upper and lower limits to the BHMF? 

• What are direct SMBH masses in QSOs and do 
they follow the established scaling relations? 

In this paper we present the observational requirements 
necessary to answers these questions. In § [2] we discuss 
the current modeling techniques that can be applied to 
large samples of galactic bulges in order to make direct 
Mm estimates. In § [3] we discuss the value of integral 
field spectroscopy in accurate M. estimates, and in § |4] 
we determine the size of telescope required to resolve 
the gravitational influence of a SMBH. § [5] discusses the 
potential sensitivities of both ground- and space-based 
facilities. In §[5] we discuss the SMBH detection abilities 
of up-coming and proposed telescopes, and §[7] sums up. 

2. DIRECT MODELS OF M, 

There are a number of direct M. estimate techniques, 
many of which are only applicable in special cases. For 
example, due to the extremely high spatial resolution 
required, proper motion studies can only be carried out 
around SagA*, and H2O Masers can only be used if the 
plane of the rotating gas is aligned very close to our line- 
of-sight. However, at present, there are two methods that 
can be applied to significant galax y populations. Models 
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Models of gas kinematics are conceptually straight- 
forward (a rotating Keplcrian disk) provided that a nu- 
clear gas disk is in-fact present and has a well defined 
inclination. In addition, the disk must be dominated by 
the gravitational influence of the SMBH, rather than by 
any inflows, outflows or turbulences. The narrow Ha 
and [Nil] emission lines are typically used to build a nu- 
clear rotation curve, however, it remains unclear as to 
whether other strong emission lines from other species, 
e.g., CIV, CIII], Mgll, [OIII], Hel, produce consistent M, 
estimates. For example, observed narrow emission lines 
are similar to the emission lines observed in planetary 
nebulae and H II regions, both of which experience strong 
non-gravitational kinematics. A final complication is in- 
troduced if there is an active galactic nucleus (AGN) 
present. AGN produce spatially unresolved broad emis- 
sion lines that must be subtracted in order to construct a 
clean rotation curve. However, the use of strong emission 
lines means that high signal to noise (S/N) observations 
can be made in a relatively short time. 

Stellar dynamical models do not suffer from many of 
the draw backs of gas kinematics (possible non-presence 
of a nuclear disk, unconstrained disk inclinations, non- 
gravitational motions) but it is also unclear as to whether 
different absorption features, e.g., Ca H & K, Mgb, CaT, 
CO band-heads, produce consistent M. estimates. The 
use of stellar absorption features means that high contin- 
uum S /N is required in order to fit the line-of-sight veloc- 
ity distributions (LOSVDs) . The low surface brightness 
of many bulges means that the required S/N can be dif- 
ficult to achieve using a practical exposure time. Finally 
stellar dynamical models are notoriously complex and 



The Future of Direct SMBH Masses 



3 



require a large amount of time to completely explore % 2 
space. In addition, there is a possible degeneracy in the 
stellar dynamical models that can produce a significa nt 
range of consistent SMBH masses (IValluri et al.ll2004h . 

Despite their respective disadvantages, gas and stel- 
lar dynamical models currently offer the only methods 
that can produce significant populations of homogeneous 
M. estimates. However, it must be noted that while al- 
most 90% of current direct estima tes have been made 
using these methods ()Grahamll2008f ). it is still unclear as 
to whether the two methods in themselves produce con- 
sistent results. Nevertheless, gas and stellar dynamical 
studies have shown three requirements of data capable 
of directly determining M. in a galactic bulges. These 
requirements are explained in more detail in the proceed- 
ing sections. First, integral field spectroscopy (IFS) must 
be used to accurately determine the stellar contribution 
to the total gravitational potential. Second, the stellar 
potential needs to be separated from the contribution of 
the SMBH by resolving the sphere of influence radius 
(rh). Finally, the S/N of the data (instrument sensitiv- 
ity) must be high enough to accurately fit line profiles 
and LOSVDs. 

3. INTEGRAL FIELD SPECTROSCOPY 

IFS is used to determine the stellar gravitational po- 
tential and provide the large number of data points es- 
sential for fully exp loring \ 2 space in dynamical models 
(jValluri et al.l l2004) . However, the value of IFS stretches 
beyond its use in dynamical models to defining the prop- 
erties of the bulge itself. A problem with the fitting 
of the M» — relation concerns the values of a * used 
(jFerrarese fc Merritt] l2000l : iGebhardt et al.ll2000fl . The 
size and shape of the spectroscopic aperture used to ob- 
serve the LOS VP has an impact on the fitted parameters 
(|Batcheldor et all 120051) . Therefore, IFS is an essential 
tool in estimating a robust M, — <r* relation as any 2D 
aperture can be used to determine cr* after the data has 
been collected. 

An integral field unit (IFU) optimized to detect 
SMBHs will be able to spatially over-sample rv In ad- 
dition, the IFU must have a field of view able to spa- 
tially sample out to an effective radius (typically ~ lkpc, 
iMarconi fc HuntJl2003h to accurately define the gravita- 
tional potential of the host bulge. The spectral pixel size 
of the IFU could vary across the field of view. Maximum 
spatial resolution is required at the center of the array, 
while spatial sampling over a larger area can be used at 
greater radii in order to collect enough S/N from lower 
surface brightness regions. Finally, the IFU must be able 
to sample both ends of the BHMF, e.g., IMBHs in globu- 
lar clusters and dwarf spheroidals (cr* » 10km s _1 ), and 
HMBHs in QSOs. Therefore, it must have high spectral 
resolution (5ft = X/5X ~ 30,000) and be able to detect 
the stellar continuum of QSO hosts while not saturating 
the detector at the QSO nucleus. 

Emissions from QSOs dominate those from the host 
galaxy. However, an IFU could directly exclude contribu- 
tions from active nuclei through the use of a micro mirror 
array. Alternatively, the active nucleus and host could 
simply be observed usin g an extreme d ynamic range ar- 
ray such as the Reticon ()Cizdziellll990l ). These detectors 
have full wells capable of holding 10 9 electrons before af- 
fecting adjacent pixels. Charge injection devices, such as 
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Fig. 1. — Angular diameter distance (Da) as a function of red- 
shift assuming the standard cosmology (Q m = 0.27, Q\ = 0.73 and 
H = 73km s" 1 Mpc" 1 ). 

SpectraCam ([Bhaskara n et al.l 120 05). may also be able 
to simultaneously observe both bright and faint sources 
in the same field of view, without compromising the in- 
tegrity of the detector, or the science image, as individual 
pixels are read out when the limit of the full well is ap- 
proached. In addition to observations of active galaxies, 
such detectors will also be valuable asset in the study of 
many other important astrophysical objects, and prop- 
erties, that require a high dynamic range, e.g., the low 
end of the stellar initial mass function, supernova ejecta, 
stellar debris disk and extra-solar planets. 

4. REQUIRED SPATIAL RESOLUTION 

In considering a collapsed obje ct at the center of a dis- 
tribution of stars, iPeeblesI (|1972| ) first derived a charac- 
teristic radius given by Equation [TJ This radius, which 
has since become known as the "sphere of influence", 
has been used to estimate the spatial scale at which the 
potential of the SMBH dominates over that of the host 
galaxy. 



GM. 



(i) 



Spatially resolving is considered to be a require- 
ment for detecting the dynamical s i gnatu r e of a SMBH 
(iFerraresd 120021 : IMarconi fc HuntJ 120031 : IValluri et all 
l2004f ). and therefore directly measuring M # . It then fol- 
lows that the M, — cr* relation can be used to determine 
the ability of a telescope to resolve Vh given its diffrac- 
tion limit (#£>)■ However, to investigate cosmic evolution 
of the SMBH scaling relations, M, must be determined 
across a large range of redshift, i.e., cosmological affects 
must be taken into consideration. 

When considering spatial resolution as a function 
of redshift, the cosmological angular diameter distance 
{Da) is of fundamental importance. It is the ratio of 
physical size to the angle subtended on the sky, i.e., the 
SMBH sphere of influence radius to the diffraction limit 
of the telescope. Therefore, in determining Da we can 
calculate the relation between the angular size of and 
redshift, and subsequently the ability of a diffraction lim- 
ited telescope to make direct determinations of M,. 

The tangential co- moving distance {Dm) is simply re- 
lated to Da by Equation [2] (where z is redshift). 
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Fig. 2. — SMBH detection requirements as a function of redshift. 
Each line marks the different selected spatial resolutions, given by 
the inset, needed to resolve r^. f0 8,8 Afo SMBHs can be seen by a 
CK'005 diffraction limited telescope at all redshifts. 

In a flat Universe (f^ = 0.00) Dm is equal to Dg (the 
radial co- moving distance). Therefore, following iPeebles! 
(1993) it can be shown that Da is given by Equation [3] 



D A = 



H (l + z) 



(n M (i + z') 3 + nxy^dz' (3) 



Assuming the standard cosmological model, with 
M = 0.27, fix = 0.73 and H = 73km s" 1 Mpc" 1 , we 
can see in Figure Q] how Da varies as a function of red- 
shift. This demonstrates a turnover in Da at a redshift 
of 1.6 (D A = 1700 Mpc rad^ 1 ). Therefore, if an object 
can be spatially resolved at this turnover redshift, it will 
be resolved at all higher redshifts. 

The M, — cr» relation can be rewritten as: 



M. 



10 a | — 



(4) 



where a = 200km s" 1 (jFerrarese fc Fordll2005l ). This 
can be rearranged for er 2 and substituted into Equation|4] 
to show: 



D A 9 D = GM, 



10 Q 



2/0 



(5) 



where D aOd = fh- Consequently, we can solve for M. 
and derive Equation^ Therefore, for a given Od, we can 
calculate the values of M. that can be directly deter- 
mined as a function of redshift. Furthermore, assuming 
a = 8.22 and (3 = 4.86 (jFerrarese fc Fordl [2005). Figure|2] 
can be produced to demonstrate the resolution require- 
ments of a direct SMBH mass estimator. It can be seen 
that e D < O'.'Ol is required to resolve M. = 1O 9 3 M , and 
9 D < 0'.'005 to resolve M. = 1O 8 8 M , at all redshifts. 
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~D A D 


( 4 )} 


G 





-1(1-2/0)- 



(G) 



Given a specific spectral feature at a particular rest 
wavelength, Figure [5] can be simply related the physi- 
cal diameter of a telescope capable of observing M, (at 
a given redshift) by using the diffraction limit equation. 
Figure[3]demonstrates these relations considering the Ha 
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Fig. 3. — Required telescope diameters to achieve specific resolu- 
tions as a function of redshift. [Upper panel] The redshift evolution 
of Ha (6563A rest frame). [Lower panel] The redshift evolution of 
Lya (f2f6A rest frame). 

and Lya emission lines that could be used for gas kine- 
matical M. estimates. These two lines are chosen to 
demonstrate the impact of being able to observe features 
in the UV verses the optical, i.e., ground-based verses 
space-based (see §0- 

Figures and [3] can be used in tandem to determine 
what size telescope is required to observe a specific mass 
SMBH at a specific redshift using Ha or Lya. From 
Figure [2] we can see that at z = 1, a resolution of 0'.'02 
is required to detect a 1O 9 8 M HMBH. From Figured 
observing this system using Ha requires a 17 m diffrac- 
tion limited telescope. However, if one were to switch 
to observe Lya, then the same observations can be made 
using a 3.1 m primary. Alternatively, the same 17 m tele- 
scope would be able to detect SMBHs down to masses of 
1O 8 ' 8 M out past z = 1.6 using Lya. Tabled] uses this 
approach to present the telescope diameters required to 
observe a range of M, across a range in redshift, using 
several different spectral features. 

5. REQUIRED SENSITIVITIES 

The same cosmology that allows the angular diameter 
distance to turn over at z — 1.6 also affects the observed 
surface brightnesses; the same flux is now distributed 
over a larger solid angle. The luminosity distance [Dl] 
is related to the angular diameter distance by Equation[7] 
(Hogg IT999h . Subsequently, surface brightness drops off 
rapidly in a standard cosmology. 



D L = D A {l + zf 



(7) 



To assess how this cosmological dimming will affect M, 
estimates we take M87 as an example. As M87 is a giant 
elliptical, it represents the class of galaxy that is expected 
to host t he most massive SMBHs , i.e., brightest cluster 
galaxies ()Dalla Bonta et all 2009). Therefore, based on 
resolving r^, such giant ellipticals offer the chance to di- 
rectly determine the mass of the highest redshift SMBHs. 
In addition, M87 type objects provide a challenge to both 
stellar dynamical models, due to low surface brightness, 
and gas dynamical models, due to weak nuclear emission 
lines. 

We have estimated a 6500A continuum surface bright- 
ness of 16.0mag arcsec -2 based on the ACS F475W 
and F850LP surface brightness profiles presented by 
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TABLE 1 

Required Telescope Diameters for SMBH Mass Estimates 



M, 

(Me) 


0.01 


Lya 
0.1 1 


10 


0.01 


0.1 


Ha 

1 


Redshift 
10 0.01 


0.1 


CaT 

1 


10 


0.01 


0.1 


Hel 

1 


10 


10 5 


27 


260 


2100 


6100 


140 


1400 


11000 


33000 


188 


1800 


15000 


43000 


240 


2300 


19000 


54000 


10 6 


6.9 


68 


540 


1600 


37 


370 


2900 


8500 


48 


480 


3800 


11000 


61 


600 


4800 


14000 


10 7 


1.8 


17 


140 


400 


9.6 


94 


760 


2200 


13 


122 


980 


2900 


16 


155 


1200 


3600 


10 8 


0.46 


4.5 


36 


100 


2.5 


24 


200 


570 


3.2 


32 


250 


740 


4.1 


10 


320 


940 


10 9 


0.12 


1.2 


9.3 


27 


0.64 


6.3 


50 


150 


0.83 


8.2 


66 


190 


1.1 


10 


83 


240 


10 10 


0.030 


0.30 


2.1 


7.0 


0.16 


1.6 


13 


38 


0.21 


2.1 


17 


49 


0.27 


2.7 


21 


62 



Note. — Telescope diameters (in meters) required to resolve specific black hole masses at specific redshifts, using specific 
spectral features. This assumes the telescopes are producing diffraction limited observations at all wavelengths, and does not take 
into account atmospheric absorption, i.e., a ground-based telescope cannot detect Lya until z>1.63. 
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Fig. 4. — Cosmological dimming of an M87 type object. Colored 
lines mark the sky backgrounds at 1.3/xm and the effective sensitivi- 
ties of a 30 m ground-based telescope, a 16 m space-based telescope, 
and JWST+NIRSpec assuming a 250 ksec exposure time. 



iFerrarese et al.l ((2006b). In addition, we have taken the 
Ha line flux of 0.8 x 10~ 15 e rgs s -1 cm" 2 as de termined 
from STIS observations by iSabra et alj (|2003l ). These 
estimates were converted to absolut e values assuming a 
distance modulus to Virgo of 31.09 (|Tonrv et al.| [2001). 
The apparent magnitudes were then rescaled, assuming 
a constant vL Vl with luminosity distance as a function 
of redshift. The relation of these apparent magnitudes 
as presented, as a function of redshift, in Figure [4] 

In § [5] we compare the abilities of both ground- and 
space-based telescopes to make direct SMBH mass es- 
timates. Therefore, for comparison, in addition to the 
continuum and line fluxes for M87 as a function of 
redshift, we over plot the sky backgrounds as expe- 
rienced at Mauna Kea 1 using the solid red line and 
at the L2 Lagrangian point using the solid green line 
([Windhorst et al. 1 12001). In addition, we plot the theo- 
retical 1.3/xm sensitivities of a diffraction limited 30 m 
ground-based telescope (dashed red line), a 16 m space- 
based telescope (dashed green line) and JWST+NIRspec 
(dashed blue line). These theoretical sensitivity limits 
have been calculated based on achieving S/N = 10 of a 
C//15 extended source in a 250 ksec exposure (M. Post- 
man private communication). The spectral resolution is 

= 2000, except for JWST where NIRSpec (5ft = 2700) 
is used. In all cases, the 16 m space-based telescope out- 



performs the 30 m ground-based facility, especially in the 
sky background levels. 

6. THE FUTURE OF M. ESTIMATES 

In § [4] and § [5] the abilities of a ground-based 30 m 
telescope, and a space-based 16 m telescope were briefly 
compared. The 30 m telescope represents the next gen- 
eration of extremely large ground-based facilities such as 
the Thirty Meter Telescope (TMT), the 24.5 m Giant 
Magellan Telescope (GMT) and the 42 m European Ex- 
tremely Large Telescope (E-ELT), that might be used 
to make the next step forward in our understanding of 
M m . Each of these will provide the extremely high sensi- 
tivities required by stellar dynamical models. The 16 m 
telescope represents the possible future of UVOIR space- 
based observatories, the size of wh ich are governed by 
future payload launch abilities, e.g.. iPostmanl (|2009l ). 

There are some other notable upcoming facilities 
that will potentially provide further interesting data on 
SMBH mas ses. The Lase r Interferometer Space Antenna 
(LISA, e.g..lHughes 2003]) will be sensitive to the gravita- 
tional wave signatures of in-falling, and coalescing binary 
SMBHs. The parameters of the system (e.g., total mass) 
can be theoretically recovered from the periodic space- 
time strain. However, even with clearly identified signa- 
tures of binary SMBHs LISA will not be able to give de- 
tails on the surrounding host galaxy, nor accurate enough 
co-ordinates for follow-up studies with other facilities. 
Another interesting possibility is to use the Atacama 
Large Millimeter/Sub- millimeter Array (ALMA) to map 
the 2.6 mm CO kinematics. In its largest base-line con- 
figuration of 1 4.5 km, ALMA will ha ve a spatial resolu- 
tion of 0"036 (|Peck fc Beaslevl 12008ft . the equivalent of 
a 4.2 m diffraction limited telescope at 600 nm. The fi- 
nal upcoming facility that may provide a step forward 
in Mm determinations is the James Webb Space Tele- 
scope (JWST). In combination with NIRSpec, the 6.5 m 
aperture will be able to provide diffraction limited inte- 
gral field and long-slit spectra (5ft ~ 3000) from 1-5/zm. 
Therefore, Paa (1.88/im) and the CO band-heads (1.5- 
4.7/im) will be observable for possible gas and stellar dy- 
namical models. Unfortunately, at these wavelengths the 
diffraction limit of JWST offers no spatial advantage over 
observations with HST. However, for 15 mag arcsec -2 ex- 



tended sources at 8561A, STIS requires > 540 minutes to 
gain a single S/N~ 50 spectrum. Assuming NIRSpec has 

1 e.g., http: //www2.gemim.edu/ sciops/telescop es-and-sites/observin g-aoaiiMd)arceffiekint^ r to STIS, the same observation could 
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Fig. 5. — Composite QSO spectrum at several redshifts as compared to the atmospheric relative transmission from Hawaii (grey shaded 
area). At z=0.07 He* becomes hampered by strong H2O water absorption at 0.7/xm. At z=0.52 Ha will be observable using NIRSpec on 
JWST at 1.0/im. At z=1.62 Lya becomes visible from the ground. 

be made in 80 minutes, making JWST significantly more 
efficient for absorption line spectroscopy. This gain in ef- 
ficiency will not translate to the emission lines, however, 
due to the relative line strength of Ha to Paa. 

Despite the possible advances offered by ALMA and 
JWST, it must be noted that there have yet to be any 
significant attempts to confirm that SMBH masses de- 
rived from gas or stellar dynamics, or from multiple emis- 
sion lines (e.g., Lya, Ha, Paa, [OIII]) and absorption 
features (e.g., Mgb, CaT), produce consistent results. 
Different spectral features will be affected by different 
issues in different ways. For example, in the case of 
absorption lines, AGN continua will fill shorter wave- 
length features, careful sky subtraction will need to be 
performed for longer wavelength features, and different 
features contend with template mismatching with differ- 
ing levels of success. There are a few cases, however, in 
which independent, i.e., from different authors, SMBH 
mass estimates have been made from both stars and gas, 
or by using different spectral features (Tab. \2§- Unfor- 
tunately, no significant conclusions can be drawn from 
such a small sample. There are a number of additional 
cases in which authors have attempted to reconcile gas 
and star estimates d irectly. For example, in NGC 3379 
iShapiro et al.l (|2006f ) find gas and stellar dynamical esti- 
mates to be in agreement (although strong non-circular 
motions are detected in the gas disk.), and in NGC 4335 
IVerdoes Kleijn et al.l (|2002h find it difficult to match gas 
and stellar dyn amical mass estimates. To add to the 
confus ion, both iKormendv et alj (fl996) Ems ellem et al.l 
(1999) both use stellar models (using the same spectral 
features) to derived M, — 10 93 and M, = 10 90 , respec- 
tively, in NGC 3115. 

While it may be possible to mitigate the effects of many 
of these issues, the most significant obstacle that must 
be faced by all ground-based observations remains the 
atmosphere. The Earth's thin protective blanket influ- 
ences spatial resolution, limits spectral coverage, and af- 
fects the required sensitivity of ground-based instruments 
with transmissions that vary with wavelength. In addi- 
tion, red-ward of ~ 0.7/im the sky begins to significantly 
interfere with ground-ba sed observations due to strong 
OH emission lines (e.g.. [Osterbrock et al.lll997f ). The 
removal of these sky lines can introduce significant over- 
heads to observing programs. For example, the Nod and 



Shuffle technique (jGlazebrook k Bland-HawthornlfeOOlh 
removes sky lines with high precision, at the expense of 
doubling the effective exposure time. The only way to 
over come these issues is to place an instrument in space, 
however, limited spatial resolution is now being success- 
fully addressed by adaptive optics (AO) systems. 

At present AO systems are limited to the near and mid- 
IR, i.e., > 1/Ltm, therefore, albeit for the significant dif- 
ferences in sky background, ground-based AO equipped 
facilities should be able to perform M, that are compa- 
rable to JWST. However, it must be noted that due to 
guide star requirements, e.g., laser guide stars and the 
availability of natural guide stars for low-order tip tilt 
corrections, full sky coverage is not currently possible for 
AO systems. In addition, further pointing restrictions 
are placed on laser guide star positioning due to possi- 
ble beam collisions between multiple systems, and due 
to Federal Aviation Administration and US Space Com- 
mand restrictions on laser projections with respects to 
flight crew distraction and satellite interference. Finally, 
AO systems must be able to produce spatially and tem- 
porally stable diffraction limited observations to enable 
deep exposures over a wide field of view. 

With careful planning and advancements in future gen- 
erations of AO systems, it may be possible to allevi- 
ate many of these issues and potentially produce diffrac- 
tion limited ground-based observations at optical wave- 
lengths. However, regardless of the future abilities of AO 
systems, ground-based observations will always be lim- 
ited by atmospheric transmission. In order to demon- 
strated these limitations, Figure [5] shows a composite 
QSO spectrum as compared to an atmospheric trans- 
mission model. The QSO spectrum combines SDSS data 
(IVanden Berk et al.ll2001l) w ith th e near-infrared (NIR) 
data of lGlikman et al.l (|2006). The lGlikman et all (|2006l ) 
data has been reduced by 0.52 flux units and chopped at 
8556A to make a smooth transition with the SDSS data. 
The emission lines that could be used by gas kinemat- 
ics models are marked, although we note that the host 
bulge absorption features that trace the stellar dynamics 
at -3950A (Ca H&K), ~5175A(Mgb), ~8570A (CaT) 
and from 1.5-4.7/mi (CO bands) will not be seen on the 
scale used. The atmospheric mod el was generated us- 
ing MODTRAN (jBerk et al.lll999h assuming conditions 
typically found at observatory altitudes on Hawaii. 
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TABLE 2 

Independent Black Hole Mass Estimates 



Name 




Mm 


Rcf. 


Mm 


Rcf. 


Consistent? 


NGC 3227 




x 1O 7 M 


s-1 


2 0+ 10 
z ' u -0.5 


x 10 7 A/ o 


g-2 


Yes 


NGC 4151 


6.5(±0.7) x 1O 7 M 


s-3 


3+1 x 


1O 7 M 


g-1 


No 


NGC 5128 


0.7 - 


- 1.1 x 10 s Mq 


g-4 a 


B±? x 


10 7 M Q 


g-5 b 


Yes 



Note. — Independent and peer reviewed direct black hole mass estimates using 
gas and stellar dynamics, or using different spectral features. References are "s" 
for stellar dyna mics, "g" for g as kin emat ics, (a) for [SHI] AA9071 , 9533 A, (b) for H2 
(2.122 itm), (1) IDavies et al.l J200l) ; (2) Iffic ks & Malkan (2008); (3) Onk en et al.l 
pOOW) : (41 IMarconi et al.l d200gT ~T5') INeumaver et al.l H2007I) 



Lya 

CIV 
CIII] 
Mgll - 
CaK 
CaH 
[OIII] 

Mgb 
Ha 

CaT 

Hel 

Pa(3 

Paa 



1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 
l|im AO Limit 



1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 
6000A AO Limit — 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 
No AO Limit 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Lookback time, t. / 1.. 

L H 

Fig. 6. — Visibility of emission and absorption features, between 800A and 3.0/jm in the rest frame, that can potentially be used to 
determine M, , as seen by any ground- based telescope, as a function of look-back time. The presence of a solid line indicates that the 
spectral feature is visible from the ground. Three adaptive optics dependent situations are shown. [Left panel] No advances in AO. [Middle 
Panel] Diffraction limited observations at 6000A. [Right panel] Diffraction limited observations all all wavelengths. 



The QSO spectrum is shown at z=0.00 and z=0.07, 
when Ha becomes hampered by strong H2O water ab- 
sorption at 0.7^m. In addition, z=0.52 is shown as Ha 
will be observable using NIRSpec on JWST at 1.0/xm. 
Paa is observable by NIRSpec at z=0 but the diffraction 
limit of JWST at 1.88^m is 0"07, the same as HST at 
Ha. Figure O also shows the QSO spectrum at z=1.63, 
when Lya starts to suffer less than 50% transmission loss 
at 3200A. 

To more precisely determine which spectral features 
are observable from the ground, and at what look-back 
time, in Figure [6] we have simply multiplied the QSO 
spectrum by the relative atmospheric transmission be- 
tween 800A and 3.0/im. If a solid line is present, then 
the indicated spectral feature is observable at the indi- 
cated look-back time. 

In comparing a 16 m to a 30 m telescope we have con- 
sidered the difference in limiting magnitudes between the 
two, according to: 



Am = 5 log 



di 



(8) 



where d\ = 30 m and di- 



16 m, and therefore Am=1.37 
magnitudes (jKitchinl [2003). The value of Am is the 
equivalent of knowing the change in magnitude as a re- 



sult of atmospheric extinction. It is simple to show that 
Am = fc(A)secZ, where fc(A) is the wavelength depen- 
dent extinction coefficient, and Z is zenith distance. In 
this case, we know that the zenith distance is zero and 
that fc(A) = 1 as it is taken into account by the rela- 
tive atmospheric transmission model. Therefore, we can 
see that the inverse of Am is equal to the relative at- 
mospheric transmission (0.73) that a ground-based 30 m 
can suffer from before it's aperture advantage is lost to a 
space-based 16 m telescope. Consequently, in Figure El 
if the atmospheric transmission is less than 0.73 then the 
spectral feature is considered blocked from ground-based 
observations. 

Figure [6] shows three different cases for the future of 
AO systems. In the first case, there is no significant ad- 
vance in AO technology, and observations are only con- 
sidered to be diffraction limited at 1.0/im. In the second 
case, AO systems are assumed to be able to produce 
diffraction limited observations down to 6000A. In the 
final case AO systems are unlimited and produce diffrac- 
tion limited observations at all wave lengths. 

Look-back time (e.g., iHogd [19991 is plotted in Fig- 
ure [6l assuming flk = 0.00, as it more directly shows the 
cosmic era over which M m can be measured. As demon- 
strated by Equation the look-back time, is similar 
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SMBH Detection Abilities 
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Fig. 7. — Required telescope diameters to achieve specific resolu- 
tions as a function of redshift. [Upper panel] The redshift evolution 
of Hel (1.1/im rest frame). [Lower panel] The redshift evolution of 
CaT (~8567A rest frame). 

to Equation [3] and proportional to the Hubble time, tn- 

dz' 



tj 



H 



o (i + z')^n M (i + z'^ + n x 



(9) 



Figure [5] is structured so that the shortest wavelength 
spectral features, and therefore the features that offer the 
greatest spatial resolution per aperture, are at the top. 
In addition, the closest, and therefore most well resolved, 
SMBHs will be found on the left hand side of each panel. , 
For reference, the Da turnover at z—1.6 (Fig. [T]) is at a 
look-back time of tL—Q.ltn- In all AO cases, the top 
left areas of each panel in Figure [5] (that can be filled by 
space-based observations and offer the greatest spatial 
resolution) remain empty. 

However, Figure [5] demonstrates that multiple spectral 
features can theoretically be used to determine M, across 
a significant cosmic era. With unlimited AO abilities, Ha 
and Mgb can be used to make gas and stellar dynamical 
models across ~ 90% of cosmic history. In addition, Ha 
can be used to cover the same era if AO systems can 
provide diffraction limited observations at 6000A. Stel- 
lar dynamical models will, with 6000A diffraction lim- 
ited observations, still be able to cover this era by using 
a combination of Mgb and CaT estimates. Using the 
current AO abilities of 1.0/^m, then ground-based obser- 
vations will need to use a combination of Hel, CaT, Mgb 
and Ha to cover a significant portion of cosmic history. 
In Figure [7] we demonstrate the size of aperture required 
to achieve these observations, for Hel and CaT, using the 
same approach as for Figure [3] 

For comparison with Figure [31 i.e., using Ha and Lya, 
Figures [3] and [7J can now be used together in order to 
determine what size telescope is required to observe a 
specific mass SMBH, at a specific redshift, using Hel or 
CaT. As seen in Figure [21 switching between spectral 
features dramatically affects SMBH detection efficiency. 
For example, to make the same observations as a 16 m, 
i.e., Ha at z=l with a spatial resolution of 0'.'02 and a 
mass of 10 9 ' 8 Mo, a 27 m diffraction limited telescope 
is needed using Hel. Switching to observe CaT then re- 
quires a 21 m primary. The same observations for Paa 
(which has a lower line flux than Hel) require a 47 m 
diffraction limited telescope. Therefore, a ground-based 
30 m telescope using existing AO abilities, will be able to 



Telescope 


Redshift ran; 




Min 


log 


Mm 




Lya 


Ha 


Hel 


Lya 


Ha 


Hel 


HST 


0.0-7.5 


0.0-0.6 




7.0 


8.6 




JWST 




0.5-6.6 


0.0-2.8 




8.5 


8.5 


16SB 


0.0-7.5 


0.0-0.6 




6.0 


7.1 




TMT(a) 


1.6-4.8 


0.0-0.1 


0.0-0.2 


8.4 


6.7 


7.1 


TMT(b) 


3.9-4.8 


0.0-0.1 


0.0-0.2 


8.8 


6.7 


7.1 



Note. — SMBH detection abilities, including redshift cov- 
erage as a function of prominent emission lines, for current and 
pending telescopes. 16SB: 16 m space-based telescope. TMT: 
30 m ground-based telescope, diffraction limited at (a) 3200A 
and (b) 6000A. "z range" : the continuous observable redshifts 
for ffi>T+STIS (1150-10,300A), and JWST+NIRSpec (1.0- 
5.0/im). We give 16SB an instrument equivalent to STIS, and 
limit TMT to 3200A - 1.3/im. "Min. logM.": the minimum 
Mm detectable at lOOMpc. 



theoretically determine the masses of high mass SMBHs 
and HMBHs out to high redshifts using Hel. However, 
as seen in Figure [5l the line strength of Hel is signif- 
icantly weaker than higher energy emission lines. Hel 
will therefore more rapidly succumb to the background 
limits shown in Figure [U 

Table [3] summarizes the SMBH detection abilities of 
several space-based telescopes and a 30 m ground-based 
telescope that can achieve diffraction limited observa- 
tions from 3200A (where Lya suffers from less than 50% 
atmospheric transmission loss) and 6000A, to 1.3/xm, 
where atmospheric transmission again impedes contin- 
uous redshift coverage. For the minimum detectable 
SMBH masses, a distance of lOOMpc is chosen in or- 
der to cover the Coma cluster, where a large range of 
SMBHs are likely to reside. Considering only the gas 
kinematics technique, the one most suited to high red- 
shift M m estimates due to emission line strengths, it can 
be seen that models using the strong hydrogen recombi- 
nation lines are best supplied data by space-based obser- 
vatories. Such facilities provide continuous redshift cov- 
erage using Lya and can therefore observe significantly 
smaller SMBHs out to significantly higher redshifts. Due 
to atmospheric transmission, at the distance of the Coma 
cluster, the 30 m ground-based telescope can only reach 
down to SMBHs of 1O 6 7 M compared to the 10 6 Af Q 
of the 16 m. 

In addition to the data in Tabled we can see that by 
using Ha the 6 D of HST can detect 10 10 M Q HMBHs out 
to a maximum distance of z=0.15. However, a sensitivity 
of ~ 25 mag arcsec" 2 is required to detect an object 
similar to M87 at this distance (Fig.g]). Therefore, HST 
can only determine M, in the brightest most massive 
QSOs at this distance, of which there are few. Using 
Lya, HSTs distance limit changes to 10 9 7 M Q SMBHs 
at z=0.6. JWST will not be able to observe Ha until 
z=0.6, but will be able to observe Hel, Pa/3 and Paa 
across a large redshift range. The longer wavelengths 
of these features, however, negates the increase in the 
JWST aperture as compared to HST, and results in no 
gain in spatial resolution. The JWST advantage will 
come from the greater collecting area. Naturally, a 16 m 
space-based mission would not be as limited as HST and 
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JWST. It would be able to detect 10 9 8 M o SMBHs to 
z = 1.0 using Ha, and 10 9 4 M Q SMBHs to z = 10 using 
Lya. It would also allow a search for IMBHs in the local 
group out to 4.5 Mpc. 

Irrespective of the SMBH detection abilities of space- 
verses ground-based telescopes, it is clear that the consis- 
tency of M, modeling techniques, using multiple spectral 
features, is vital for continued investigations into the role 
of SMBHs in galaxy formation and evolution. For exam- 
ple, for the abilities listed in Table [3] to be reached, it 
must first be shown that both Lya and Hel can be used 
to determine M.. In addition, as gas disks are not ubiq- 
uitous at the nuclei of galactic bulges, gas kinematical 
models must also be reconciled with stellar dynamical 
models using at least both Mgb and CaT. The next step 
in M, investigations must therefore be to produce consis- 
tent SMBH masses, from different methods and spectral 
features, to pave the way for more complete SMBH inves- 
tigations using both space- and ground-based facilities. 
Space is required for Lya M, estimates, but ground-based 
AO assisted observations of Hel and the CO band-heads 
will be useful, provided the aperture is large enough to 
compensate for the longer wavelengths and lower line 
fluxes. Due to the large wavelength coverage (1150- 
10,300A) #ST+STIS provides an excellent opportunity 
to complete a M, calibration for Lya and both Mgb and 
CaT. In addition, the calibration of NIR lines would be 
an excellent program for JWST+NIRspec. 

7. CONCLUSIONS 

Our understanding of the co-evolution of SMBHs and 
galactic bulges has been driven by the discovery of scaling 
relations between the two. However, there are a number 
of fundamental concerns with these relations that will 
have an impact on galaxy formation and evolution sce- 
narios. The future of direct M, determinations is there- 
fore likely to be dominated by five key questions. Are 
IMBHs and HMBHs consistent with the bulge scaling 
relations defined by SMBHs? Are bulge scaling rela- 
tions linear? Have bulge scaling relations evolved to their 
present form? Are there upper and lower limits to the 
BHMF? What are the direct M. estimates in QSOs, and 
do they follow the bulge scaling relations? It is therefore 
necessary to expand our abilities to determine the masses 
of a range of SMBHs across as wide a cosmic history as 
possible. 

We have considered the requirements for an effective 
SMBH detector that uses current gas and stellar dynam- 
ical models and techniques to derive M # . Cosmological 
models have been coupled with the sphere of influence ar- 
gument to determine the aperture sizes needed in order 
to resolve SMBHs of specific masses at a range of red- 
shifts. In addition, we have considered the sensitivities 
that such facilities will need. 

It has been demonstrated that the limits of M. esti- 
mates from HST are being reached and cannot be signif- 
icantly expanded by JWST. Therefore, additional facil- 
ities are needed to directly determine whether IMBHs, 
HMBHs and QSOs follow the same scaling relations as 
typical galactic bulges, and whether locally established 
scaling relations are linear and cosmically constant. 

A 30 m ground-based optical- NIR observatory (e.g., 



the forthcoming TMT) has been directly compared to a 
UVOIR 16 m space-based observatory. For either facil- 
ity to be an effective SMBH detector, current modeling 
techniques must be directly compared with each other to 
ensure consistent mass estimates from different models 
using different spectral features. 

If consistency can be established, then a space-based 
16 m telescope can theoretically be used to estimate 
M, using Lya out to z=10, or across 96% of cosmic 
history (depending on the instruments, and limited to 
M. > 1O 9 ' 4 M ). This facility will also be able to de- 
tect M, ~ 1O 6 O M0 at a distance of the Coma Cluster, 
and M. - lO 4 5 Af in the Local Group. A 30 m ground- 
based telescope will need to use multiple spectral features 
in order to cover the same amount of cosmic history, even 
in the face of significant advances in AO systems. Due 
to limitations imposed by atmospheric transmission, the 
30 m can theoretically be used to estimate M, using Lya 
from z=1.6 to z=4.8, or across 20% of cosmic history 
(depending on the instruments). For increased cosmic 
history coverage, longer wavelength spectral features will 
need to be used at the cost of SMBH detection efficiency 
The 30 m will also be able to detect M. ~ 10 6 7 M Q at 
a distance of the Coma Cluster, and M, ~ 10 4 2 M Q ) in 
the Local Group. 

However, the abilities of the ground-based 30 m are 
critically dependent on future advances in AO systems. 
As AO abilities increase toward the NIR, so does the 
effective spatial resolution, and therefore the SMBH de- 
tection efficiency. Considering current AO systems, the 
30 m is limited to observing Hel. In this case, the 30 m 
will be able to detect M. ~ 10 71 M Q at a distance of the 
Coma Cluster. In addition, emission lines at these wave- 
lengths are significantly weaker than features blue-ward 
of Ha. As a consequence, the advantages of the large 
aperture are lost as it must integrate for longer to gain 
the same S/N that a 16 m would have observing Lya. 
The ultimate advantage for a space-based telescope, re- 
gardless of AO systems, then becomes the limitations 
to sensitivities as determined by sky backgrounds. As 
shown in Figure |4j the magnitude limits for an object 
such as M87 are reached very quickly from the ground. 
It will then become incredibly challenging, considering 
all the potential overheads (such as nod and shuffle), to 
gain the required S/N for dynamical models in a reason- 
able amount of observing time. 
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